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The pulse tube engine represents the thermodynamic inversion of the pulse tube refrigerator used in 
cryogenic cooling applications. It has a high potential to be used as a prime mover for the conversion of 
low grade waste heat into mechanical or electrical energy. This paper describes an analytical analysis of 
the pulse tube engine based on a zero-dimensional model. During compression and expansion, the 
engine components are considered as isothermal with characteristic temperatures. At the piston's dead 
centers, a thermal relaxation model is used to switch between these temperatures. Analytical relations 
for the pV— work developed by the pulse tube engine and its efficiency are derived. The irreversible 
nature of the pulse tube engine is studied by calculating the entropy production in the components. 
Furthermore, the thermodynamic cycle is investigated analytically under variation of design features and 
operating conditions. The results are compared to prior numerical studies. The minimal temperature 
ratio above which the engine provides a work output is derived analytically and compared to experi¬ 
mental observations. Fundamental characteristics and application limitations of the pulse tube engine 
are disclosed. An upper limit for the efficiency of the pulse tube engine is derived theoretically and 
confirmed experimentally as well as through numerically calculations. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

In the 1960’s Cifford and Longsworth invented the basic pulse 
tube refrigerator [1], These coolers have been research and devel¬ 
opment objects for the last 40 years, e.g. Refs. [2,3], and are 
employed in cryogenic technologies to provide as one-stage re¬ 
frigerators temperatures around the condensation temperature of 
nitrogen (77.2 1< at atmospheric pressure) or as two-stage re¬ 
frigerators temperatures down to the condensation temperature of 
helium (4.2 K at atmospheric pressure). 

The reversing of the pulse tube process enables the design of a 
simple prime mover — the pulse tube engine — for the conversion of 
externally supplied thermal energy into mechanical energy or, us¬ 
ing an electrical generator, into electrical energy. 

The basic engine design is shown in Fig. 1. A compression/ 
expansion cylinder supplies a pressure swing and gas shifting to the 
connected system consisting of an almost adiabatic buffer space 
(pulse tube) and an isothermalizer (regenerator). Both components 
are fluidically linked via a gas heater (hhx) at the temperature 7),. 
providing the heat input. A gas cooler (chx!) at the temperature T c is 
located between cylinder and pulse tube. A heat sink (chx 2 ) at the 
temperature T c — to set the depicted material temperature distri¬ 
bution — is located at the end of the regenerator. For the sake of 
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simplification, it is assumed that the material temperature is 
spatially and temporally constant in the heat exchangers and has 
the shown temporal constant, linear distributions in the pulse tube 
and regenerator. According to the numerical and experimental re¬ 
sults of [4], the linear material temperature distribution in the 
regenerator is a good assumption, whereas the material tempera¬ 
ture distribution of the pulse tube is slightly non-linear. This non¬ 
linearity has only a second order influence on the heat transfer 
conditions in the pulse tube, and is thus neglected. 

The conversion of compression/expansion work {pV— work) into 
mechanical work and vice versa is carried out by the resulting 
pressure force on the piston within the cylinder, which is located at 
cold temperature level. For the conversion of mechanical energy 
into electrical energy, an electrical generator can be joined to the 
cylinder/piston configuration via a crank gear mechanism. 

The mode of operation of the pulse tube engine is based on a 
simultaneous compression and moving of the working gas into an 
area of high material temperature followed by a simultaneous 
expansion and moving of the working gas into an area of low 
material temperature. Therefore, the piston moving from its bot¬ 
tom dead center (BDC) to its top dead center (TDC) compresses the 
working gas and pushes gas from the cylinder via chx! into the cold 
side of the pulse tube. Passing chxi, the gas is cooled down toT c . The 
situation is illustrated in Fig. 1 by the dotted gas temperature dis¬ 
tribution in the pulse tube during the compression stroke. Mean¬ 
while, in the pulse tube, gas is compressed and shifted from the 
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List of symbols 

Sirr.reg 

entropy production regenerator 



T 

temperature 

H 

enthalpy flux 

1° 

temperature of reference state 

m 

mass flux 

T c 

heat sink temperature 

^chx 

mass flux via chx 

T h 

heat source temperature 

^hhx 

mass flux via hhx 

Tchx 

gas temperature chx 

Qhx 

gas-material heat flux 

Tcyl 

gas temperature cylinder 

s 

entropy flux 

Thhx 

gas temperature hhx 

^hx 

gas-material entropy flux 

^m.chx, 

material temperature chx! 

w 

compression/expansion (pV— ) power 

^rn,chx 2 

material temperature chx 2 

C| 

fitting parameter 

Tm.cyl 

material temperature cylinder 

Q 

compression constant 

Tm.hhx 

material temperature hhx 

C e 

expansion constant 

7m.pt 

material temperature pulse tube 

Qip 

high pressure constant 

Tm.reg 

material temperature regenerator 

C.r 

low pressure constant 

Tpt.c 

compression gas temperature pulse tube 

/ 

operating frequency 

Tpt.e 

expansion gas temperature pulse tube 

7-pt 

pulse tube length 

V 

gas temperature pulse tube 

L re g 

regenerator length 

Treg 

gas temperature regenerator 

m 

gas mass 

U 

internal energy gas 

p° 

pressure of reference state 

V° 

volume of reference state 

Pc 

compression pressure 

Vd 

displacement volume 

Pe 

expansion pressure 

Vhx 

total heat exchanger gas volume 

Pf 

filling pressure 

Wnax 

maximum engine gas volume 

P S 

gross power 

Miiin 

minimum engine gas volume 

Pn 

net power 

Vpt 

pulse tube gas volume 

A.int 

internal power loss 

kreg 

regenerator gas volume 

Pmax 

maximum pressure 

W c 

compression work 

Pmin 

minimum pressure 

W e 

expansion work 

Qchx 

cycle heat chx 

W pV 

pV— work 

Qcyl 

cycle heat cylinder 

Xpt 

pulse tube position 

Qhhx 

cycle heat hhx 

X Teg 

regenerator position 

Qin 

supplied heat 

F 

fitting parameter 

Opt 

cycle heat pulse tube 

£ 

compression/expansion factor 

Qreg 

cycle heat regenerator 

e c 

compression factor 

r c 

compression ratio 

£ e 

expansion factor 

hi 

distribution ratio 

Vth 

thermal efficiency 

s 

gas entropy 

X 

boundary condition factor 

•^irr.chx 

entropy production chx 

n 

pressure ratio 

•^irr.cyl 

entropy production cylinder 

a 

sign function 

Sjrr.hhx 

entropy production hhx 

T 

temperature ratio 

•^irr.pt 

entropy production pulse tube 

^opt 

optimal temperature ratio 


cold to the hot side of the tube. Gas leaving the pulse tube at its hot 
side via hhx takes on heat and enters the hot side of the regenerator 
with the temperature Th. The regenerator material, having the 
depicted material temperature distribution of Fig. 1, absorbs heat 
from the working gas. Thus, the heat taken from hhx is entirely 
given off to the regenerator material. The thermal capacity of the 
regenerator material is high enough to experience no significant 
temperature swing. Since cold gas was shifted during compression 
from the cold to the hot side of the pulse tube, it is able to take on 


heat from the tube’s wall material while the piston is moving 
around TDC (high pressure thermal relaxation). The pressure in¬ 
creases slightly. Now the piston moves from TDC back to BDC which 
results in an expansion of the working gas flowing from the 
regenerator via hhx into the hot side of the pulse tube. The situation 
is illustrated in Fig. 1 by the dotted gas temperature distribution in 
the pulse tube during the expansion stroke. The regenerator ma¬ 
terial gives back the heat stored during the compression phase. 
Since the gas entering hhx is already at Th, it receives only an 



Fig. 1 . Basic design of the pulse tube engine, chx! - cold heat exchanger 1, chx 2 - cold heat exchanger 2, hhx — hot heat exchanger. For the sake of simplification, the material 
temperature distribution is assumed to be linear in the pulse tube and the regenerator as well as spatial constant in the other components. The ideally assumed gas temperature 
distribution in the pulse tube during compression and expansion is illustrated using dotted lines. 
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insignificant amount of heat. Meanwhile, in the pulse tube, gas is 
expanded and shifted from the hot to the cold side of the tube. Gas 
leaving the pulse tube at its cold side via chxi is cooled and enters 
the cylinder with the temperature T c . In the cylinder the gas is 
further expanded and its temperature goes below T c . Since hot gas 
was shifted during expansion from the hot to the cold side of the 
pulse tube, it is able to give off heat to the tube’s wall material while 
the piston is moving around BDC (low pressure thermal relaxation). 
The pressure decreases slightly till the initial state is reached again. 
Absorbing heat at high pressure and rejecting heat at low pressure 
results in a higher mean pressure during expansion than during 
compression and in a release of pV— work. 

First concepts of prime movers similar to the pulse tube engine 
were presented in 1987 by Chen and West 5]. In 1995, Tailer was 
granted a US patent [6] for a heat engine based on a phase delay 
between pressure fluctuations and heat transfer (thermal lag [7,8]). 
A patent describing a device resembling the pulse tube engine was 
published in 2000 by Hofmann [9], The concept of the pulse tube 
engine was first presented in its entirety by Hamaguchi et al. [10] in 
2005, and later by Organ [11] in 2007. A number of experimental 
investigations to characterize this type of prime mover were also 
published by Hamaguchi et al. [10,12], 

According to work flux measurements of Yoshida et al. [13] in 
2009, the pulse tube engine can be classified into the thermoa¬ 
coustic [14] standing wave engine group [15—18] with an under¬ 
lying intrinsically irreversible thermodynamic cycle [19], although 
the pulse tube engine is — by comparing its working principle with 
the thermoacoustic engines investigated in Refs. [20—22] — anal¬ 
ogous to the Stirling engine, not a thermoacoustic device. 

In contrast to the Stirling engine, the pulse tube engine has no 
displacer for the provision of the phase delay between pressure 
fluctuation and volume flux, which is operationally necessary. In 
similarity with the pulse tube refrigerator, the gas dynamics and 
heat transfer conditions in the pulse tube and the regenerator 
enable the generation of this phase delay. Thus, the pulse tube 
engine operates with only one piston. Due to the difficulty in 
adjusting the mentioned phase delay by the correct assembly of 
pulse tube and regenerator, the design of the pulse tube engine is 
bound to a deep understanding of the physical phenomena related 
to its operation. In order to meet this challenge, a comprehensive 
thermodynamic analysis using a numerical model of the pulse tube 
engine was performed by Moldenhauer et al. in 2011 and published 
in 2013 [23]. There, it was shown that a broken thermodynamic 
symmetry between the pulse tube and the regenerator is the un¬ 
derlying working principle of the pulse tube engine. Thermody¬ 
namic asymmetry results from the fluidic linkage of the almost 
adiabatic pulse tube and the regenerator. Thus — in contrast to the 
Stirling engine — the regenerator is essential for the operation of 
the pulse tube engine. The regenerator of the pulse tube engine can 
be regarded as an isothermalizer, which together with the almost 
adiabatic pulse tube creates a thermodynamic asymmetry. 

An extended version of the numerical model of [23] was pre¬ 
sented by Moldenhauer et al. in 2013 [24]. The numerical results are 
compared to experimental performance data of a pulse tube en¬ 
gine, built using pressurized Helium as working fluid. The numer¬ 
ical model presented in Ref. [24] can be used to design a working 
pulse tube engine and predict its performance. 

Although the pulse tube engine is thermodynamically similar 
to the pulse tube refrigerator, its theory is still much less un¬ 
derstood than for its cryogenic counterpart. To the best of the 
author’s knowledge, no closed solution of the pulse tube engine’s 
governing equations has been obtained yet. In order to bridge 
this gap, the present work represents a contribution to the the¬ 
ory of the pulse tube engine by the development of an analytical 
model for a closed solution of the underlying differential 


equations. The aim of this research is to determine the basic 
physical characteristics of the pulse tube engine, to find analyt¬ 
ical relations for the developed pV —work and its efficiency as 
well as performance limitations. The irreversible nature of the 
pulse tube engine is studied by calculating the entropy produc¬ 
tion in the components. Furthermore, the thermodynamic cycle 
is investigated under variation of design features and operating 
conditions. Finally, the analytical results are compared to nu¬ 
merical and experimental results obtained from the numerical 
models given in Refs. [23,24] and the working pulse tube engine 
presented in Refs. [4,24], 


2. Governing equations 


The system of differential equations underlying the pulse tube 
engine is based on a control volume analysis [25,26]. Each compo¬ 
nent of the engine is regarded as one control volume with one 
representative value of each physical property depending only on 
time. The balances for mass, energy, and entropy of the working gas 
lead to the governing equations for the time-dependent physical 
properties of the working gas. The change in the internal energy of 
the material is disregarded. Since pulse tube engines are typically 
operating at low frequencies, a momentum balance is not performed. 

Fig. 2 shows the schematic representation of the control volume 
V. The internal energy of the gas is represented by U. The gas mass 
content of the control volume is m and its entropy S. The material 
has the constant temperature T m . Heat can be exchanged between 
the engine material and the working gas via the heat flux Q hx . The 
entropy flux due to heat transfer between gas and material is 
represented by S hx . The working gas in the cylinder is charged with 
the compression/expansion power (pV— power) W. 

Between the control volumes, mass fluxes m, enthalpy fluxes H, 
and entropy fluxes S are exchanged. Fluxes at the control volume’s 
left side are marked with ' and fluxes at the control volume’s right 
side with ", Heat and entropy fluxes due to heat conduction are not 
considered. All entering fluxes are given a positive sign. 

The resulting system of differential equations for the pressure p, 
the gas temperature T, and the entropy production rate dSj rr /dt in 
each engine component is given by 


dp 

dt 


(m' + m") 


p dV pdT 
V dt + Tdf ’ 


(1) 


c v pv dr 

R S T dt 


-c v T(m'+ m") + Qhx - P 


dV 

df 


+ m'cpT' + m"c p T" 


( 2 ) 



Fig. 2. Schematic representation of the control volume with the internal energy of the 
gas U, the gas mass m, and the gas entropy S. Heat and entropy exchanged between the 
gas and the material are represented by Qhx S hx . Mass, enthalpy, and entropy 
fluxes between the control volumes are denoted by m, H, and S. Left and right side 
fluxes are marked with ' and The compression/expansion power is W and the ma¬ 
terial temperature is T m . 
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and 

^L=(m' + m") (c v (in^ - l) + 

+ ri 1 '( Cp (y“ ln |o) +Rsln ^) (3) 

+ m"( Cp (L_ ln I_) +R s \n^j, 

where T°, p°, and V° are the physical properties of the gas at the 
reference state. 


heat exchanger has the uniform material temperature T m .hhx — 7h. 
Pulse tube and regenerator have the space-dependent fixed linear 
material temperature distributions 

7m.pt (xpt) = 7 C + j?-{T h — T c ) 

Ept 

and 

7 m,reg(Xreg) = 7 ), — J —~( 7 h ~ 7 c), 

L reg 


3. Parameter definitions 


A generalized characterization of the pulse tube engine process 
requires the definition of dimensionless operating and design 
parameters. 

The operating parameters are the temperature ratio x, defined 
as the ratio between the hot temperature Th and the cold temper¬ 
ature T c 


as well as the pressure ratio 77, representing the ratio between 
maximum pressure p max and minimum pressure p m j n 


77 ._ Pmax ( 5 ) 

Pmin 

Design parameters are the distribution ratio r c \, calculated as 
the ratio between the regenerator gas volume V reg and the pulse 
tube gas volume Vp t 


r d : = 


Vreg 

Vpt 


( 6 ) 


as well as the compression ratio r c , representing the ratio between 
the maximum engine gas volume Vmax and the minimum engine 
gas volume i/mm minus one 


Vmax .j _ Vpt + V re g + Vjj x + V j 
Vmin Vpt + V re g + V hx 

Vd 

Vpt + Vreg + Vhx 


(7) 


with x pt denoting the pulse tube position starting at its cold side 
with Xpt = 0 and x re g denoting the regenerator position starting at 
its hot side with x reg = 0. The length of the pulse tube is L pt and of 
the regenerator I reg . 

The gas volumes of the heat exchangers are neglected. Thus, the 
cold heat exchanger chx 2 has not been considered in this analysis. 
For pulse tube and regenerator, the (spatial) mean material tem¬ 
peratures are given by 


< 7m,pt > < 7 m ,reg > — ^ (7h + Tc) — ( T + 7 )■ 

In the isothermal analysis of the pulse tube engine the gas in the 
cylinder, the heat exchangers, and the regenerator is always in 
thermal equilibrium with the material [28]. The constant gas 
temperatures in these components are T cy i = T g . cy i = T c , 
7chx = Tg.chx, = 7c, Thhx = 7 gi hhx — Th, and 7 reg = T gireg = 
<7 m , reg >=T c (T + 1 )/2. 

Due to the cold gas being pushed into the cold side of the pulse 
tube during the compression stroke, the pulse tube’s mean gas 
temperature is below its wall material temperature while the 
compression takes place. During the expansion stroke, hot gas is 
pushed into the hot side of the pulse tube. Thus its mean gas 
temperature is above the wall material temperature while the 
expansion takes place. The situation is illustrated in Fig. 3. For 
compression and expansion, it is assumed that the gas in the pulse 
tube has a constant temperature slightly below the mean wall 
material temperature (compression stroke) respectively slightly 
above the mean wall material temperature (expansion stroke). 

Using this assumption, the pulse tube’s gas temperature is set to 

Tpt.c = (7 — £c) < 7m.pt > = (7 — £c) ^ ( T 7“ 7) (8) 

during the compression stroke and 


The displacement volume of the cylinder and the sum of the 
heat exchanger gas volumes are denoted with Vd and Vh x , 
respectively. 

4. Analytical model of the pulse tube engine 

Pulse tube engines typically operate at low frequencies of only a 
few Hz. Such slow thermodynamic processes can be approached 
using an isothermal model [27], An isothermal analysis (dT/df = 0) 
of the thermodynamic cycle decouples the system of differential 
equations (1), (2), and (3) and enables an analytical integration of 
the decoupled differential equations. 

4.1. Isothermal approach 

As already shown in Fig. 1, the cylinder has the uniform material 
temperature T mcy i = T c , the two cold heat exchangers have the 
uniform material temperatures T m chX] = 7 m c hx 2 = T c , and the hot 


7pt,e = (7 + Be) < 7 m .p t > = (7 + e e )^(t + 7) (9) 



Fig. 3. Gas and material temperature of the pulse tube during compression and 
expansion. In the isothermal model, the gas temperature in the pulse tube during the 
compression T pt c is assumed to be constant and slightly below the material temper¬ 
ature 7^, whereas the gas temperature during the expansion r pt ,e is assumed to be 
constant and slightly above the material temperature. 
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during the expansion stroke. The parameters e c and e e are in the 
range of 0 < e c < 1 and 0 < e e < 1, respectively. 

At TDC, the pulse tube’s gas temperature undergoes a thermal 
relaxation [11] from r pti c to T p t, e , and at BDC, a thermal relaxation 
from Tpt,e to 7pt,c Both relaxation processes occur due to shifting 
and mixing of the gas in the pulse accompanied by entropy pro¬ 
duction [29]. Since the relaxations take place at TDC and BDC, the 
gas volume of the entire engine is constant and the differential 
equations can be solved analytically [30]. 

Summarizing, the thermodynamic cycle is split into four phases. 
According to Fig. 4, an isothermal compression at T p t, c occurs be¬ 
tween state 1 and state 2. A high pressure thermal relaxation from 
Vc to 7pt.e takes place between state 2 and state 3 at constant 
minimum engine gas volume. An isothermal expansion at lpt,e 
occurs between state 3 and state 4 and a low pressure thermal 
relaxation from T pt , e to T p t, c takes place between state 4 and state 1 
at constant maximum engine gas volume. 

4.2. Boundary conditions 

Between the engine components, boundary conditions for all 
transferred quantities need to be assigned. 

The gas flowing between the components carries thermal en¬ 
ergy and entropy. Both are related to the mass flux and the gas 
temperature. In Fig. 5, the temperatures of the gas crossing the 
control volume borders are designated. Gas leaving and entering 
the cylinder at its right side, respectively leaving and entering the 
cold heat exchanger at its left side has the temperature 
T"yi = T'chx = T c- Gas leaving and entering the hot heat exchanger 
at its right side, respectively leaving and entering the regenerator at 
its left side has the temperature T[ 1 ' hx = T' eg = T h . Gas leaving the 
cold heat exchanger at its right side, respectively entering the pulse 
tube at its left side has the temperature T" hx = Tp t = T c . Gas 
leaving the hot heat exchanger at its left side, respectively entering 
the pulse tube at its right side has the temperature 
T" t = T^ hx = T h . During the compression stroke (1 —> 2), gas is 
leaving the pulse tube at its right side and entering the hot heat 
exchanger at its left side with the temperature T^ t = T^ = T" t c , 
whereas during the high pressure thermal relaxation (2 —> 3), gas is 
leaving the pulse tube at its right side and entering the hot heat 
exchanger at its left side with the temperature 

T'' t = T^ hx = Tp tc .T" te . During the expansion stroke (3 -> 4), 

gas is leaving the pulse tube at its left side and entering the cold 
heat exchanger at its right side with the temperature 
T" hx = Tp t = Tp te . During the low pressure thermal relaxation 
(4 —► 1), gas is entering the pulse tube at both sides. 

Since all components except the pulse tube are isothermal 
throughout the entire cycle, the denotations ' and " are 
only necessary to identify the boundary values at the pulse 



Fig. 4. pV-diagram of the idealized pulse tube engine cycle. TDC: piston's top dead 
center, BDC: piston’s bottom dead center. 



Fig. 5. Boundary gas temperatures of the engine components during the four phases of 
the cycle. T^,: cylinder, T'^, 7"^: cold heat exchanger, T^ t , T" : pulse tube, r hhx , r"^: 
hot heat exchanger, TJ eg : regenerator. The arrows between the components symbolize 
the direction of the respective mass fluxes. The mass flux via the cold heat exchanger 
m dlx is zero during the high pressure thermal relaxation. 

tube’s cold ' and hot " side. The respective pulse tube boundary 
temperatures T' and T" are assumed to be 

V = )! < Tm.pt > +(1 - A')T C 

at its cold, and 

T" = A" < T m , pt > +(1 - A") T h 

at its hot side. 

During the compression stroke, the pulse tube’s hot side exit 
temperature is T" tc = (1 - e c )T" and during the expansion stroke, 
the pulse tube’s cold side exit temperature is Tp t e = (1 + e e )T'.The 
hot side exit temperature during the high pressure thermal relax¬ 
ation is Tp tc ...Tp te = (1 - ec )T"...(l +e e )T". 

Gas is entering the pulse tube’s cold side with the temperature 
T c during the compression stroke and the pulse tube’s hot side with 
the temperature Th during the expansion stroke. Hence, the pa¬ 
rameters A' and A" can be deduced from the conditions Tp t c = T c 
and Tp t e = T h to 

y _ 2e c 

" (l-*c)(T-l) 

and 

yi _ 2 Tec 

~ (l+ ee )(T-l)‘ 

It follows 

^.e = and T" tc = l^|T c T. (10) 

4.3. Calculation of the time-dependent pressure 

Eq. (1) is used to calculate the pressure throughout the ther¬ 
modynamic cycle of the pulse tube engine. For all considerations, 
the mass fluxes via the cold and hot heat exchanger are denoted 
by m chx and rn h hx- Not considering the heat exchanger gas volumes 
leads to m' chx + m" hx = 0 and m' hhx + rii[( hx = 0, and thus 
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Khx = -m'c'hx = m C hx and m' hhx = -m'h hx = m hhx , respectively. 
Left side mass fluxes m' chx and m' hhx are positive, and right side 
mass fluxes m" hx and m dhx are negative. 

With the definitions ( 6 ) and (7), Eq. (1) becomes 


Pc(t) = Pi 


V d + C c 

^cylW + Cc 


which reaches in state 2 


dp _ R S T C . p dY/ cyl 
dt “ V cyl mchx V cyl dt 


( 11 ) 


P 2 = Pi 


Vd + Cc 

Q 


( 21 ) 


( 22 ) 


for the cylinder, 

cit = ^^ rc ( rd + 1 )( rilchx_ri,hhx) (12) 

for the pulse tube during the isothermal compression and the 
isothermal expansion, 


dp R s T pt , .... . , p dipt 

dt = -vT^ + 1)(mchx “ mhhx) + iptlF 

for the pulse tube during the thermal relaxations, and 


dp 

dt 


RsTreg r c (r d + 1) 
“Vd r A 


d^hhx 


(13) 


(14) 


for the regenerator. 

The cylinder gas volume at BDC is l/ cy i = V d and at TDC V cy i = 0. 
The resulting mass fluxes can be calculated from Eqs. (11) or (12) 
and Eq. (14). The mass flux via the cold heat exchanger is during 
the isothermal compression and the isothermal expansion 


m U d / r d 1 \ dp 

chx,c/e Pst c (r d + 1) \T r eg TptJ dt’ 

during the high pressure thermal relaxation 

dtchx.hp = 0, 

and during the low pressure thermal relaxation 


(15) 


(16) 


Thermal relaxation from state 2 to state 3 

Inserting Eqs. (16) and (18) in Eq. (13), together with the 
abbreviation 



results in 


1 dp __ 1 dipt 

pdt” T pt + C hp T p 2 t dt 


(23) 


(24) 


Integrating Eq. (24) from state 2 to state 3 leads to the pressure 


P3 = P 2 


(Qipipt.c 

(Oiplpt.e 


+ 1 ^ lpt,e 
+ l) lpt,c 


(25) 


Isothermal expansion from state 3 to state 4 

Inserting Eq. (15 ) in Eq. (11) leads to 

1 dp 1 dV cy i 

p dt - Vcy] + C e dt ’ 

with the abbreviation 

c ._ Vi ( r Tc T c \ 
e ' r c (r d + l) V d T reg T pt ,e/ 


(26) 


(27) 


_ U d dp 

mchx lP “ R S T C dt ’ 


Integrating Eq. (26) with the initial conditions V cy j(t 3 ) = 1/3 = 0 and 
p(t 3 ) = p 3 = p max leads to the pressure during the expansion stroke 


Throughout all phases, the mass flux via the hot heat exchanger 
is given by 


Pe(t) = p 3 


C e 

Vcyl(t) + Ce 


^hhx — 


Vd t d dp 

Pst c (r d + 1) Treg dt 


(18) 


The pressure during the four engine phases is calculated as 
follows: 


which reaches in state 4 


P 4 = P 3 


C e 

Vd + Ce ’ 


(28) 


(29) 


Isothermal compression from state 1 to state 2 

Inserting Eq. (15) in Eq. (11) leads to 

1 dp 1 dl/ cyl 

pdt “ l/ cyl + C c dt ’ 

with the abbreviation 
c ._ Vd ( r T c T c \ 

C ' r c (r d + l) V^reg W 


(19) 


( 20 ) 


Integrating Eq. (19) with the initial conditions V'cyi(fi) = Vi = V d and 
p(tj) = pi = p min leads to the pressure during the compression 
stroke 


Thermal relaxation from state 4 to state 1 


Inserting Eqs. (17) and (18) in Eq. (13), together with the 
abbreviation 


^ r c(t d + 1) , r d 

Mp ‘ t ' y 5 

1 c 1 reg 

results in 

1 dp _ 1 dr pt 

P dt ~ Tpt + C lp T p 2 t dt ’ 


(30) 


(31) 


Integrating Eq. (31) from state 4 to state 1 leads to the pressure 
























S. Moldenhauer / Energy 62 (2013) 285-299 


291 


Pi = P 4- 


QpPpt.e 

CipT pt ,c 


+ T'pt.C 

+ 7pt,e 


(32) 


Closing the cycle through the combination of Eqs. (22), (25), 
(29), and (32) demands 


Pi 

Pi 


V d + C c ( C hp7pt,c + l) Tpr.e C e (Qp^pte + l^Tpt.c , 

(Qip7pt,e + l) 7pt,c + *~ e (C| p Tpt,c + 1 j Tpt e 

(33) 


which is fulfilled for all values of e c and e e . Without loss of gener¬ 
ality, the ratio between e c and e e is defined as 


To : =-- ( 34 ) 

With e ■■= e e results e c = Toe. The value of e can be obtained from the 
mean gas temperature in the pulse tube during compression 


Tpt,c = ^(T' t:C + T" tc ) (35) 

or during expansion 


T pt ,e = jfr't.e + T'pt.e), (36) 

respectively. Eqs. (35) and (36) lead to the same result for 


t/tq ^ 1 

T + 1 


(37) 


4.4. Calculation of the pV—work and transferred heat 

The pV—work W pV is calculated through integration of the 
isothermal compression power W c and the isothermal expansion 
power VV e according to 

2 4 o v d 

W pV = j W c dt + J Wedt = - f p c dV cyl - j p e dV cyl . (38) 

13 V d 0 

Eq. (38) results with Eqs. (21) and (28) in 

WpV = Pmin(V d + Cc)ln(g+l) - pmaxCeln^ + lj. (39) 

The heat fluxes Q cyl , Q pt , and Q reg to the gas in the respective 
components are calculated using Eq. (2), which, for the cylinder 
becomes 


0 = c v T c m chx + Q C yi - p - m chx c p T c , (40) 

for the pulse tube during the compression and expansion stroke 

0 = c v T p t(m hhx - m chx ) + Q pt + m chx c p T" - m hhx c p T', (41) 

for the pulse tube during the thermal relaxations with the defini¬ 
tions (6) and (7) 


Psrc(r d V + 1) 'dr = CvT(rilhhx “ ™chx) + + ^chxCp T" 

~ m h hxCp T', (42) 

and for the regenerator 


0 — —CvTreghlhjjx + Q re g + rhhiaxCpTh. (43) 

Using Eq. (2), the heat fluxes in the heat exchangers are calcu¬ 
lated for the cold heat exchanger with 

Qchx = m chx c p (T"-Tc), (44) 

and for the hot heat exchanger with 

Qhhx = m hhx Cp(T h - V), (45) 

respectively. 

Solving Eq. (40)—(45) for the four phases and integrating the 
heat fluxes over the entire engine cycle leads to the amounts of heat 
transferred to the working gas or removed from the working gas in 
the components per cycle according to 


Qcyl/chx/pt/hhx — T Q-cyl/chx/pt/hhxdf 


(46) 


The amount of heat transferred per cycle in the regenerator is 


Qreg 



(47) 


4.5. Calculation of the entropy production 


To understand the power and efficiency limitations of the pulse 
tube engine, the entropy production rate is calculated using Eq. (3). 
For the cylinder follows 

%^ = -m ch x(cv( ln£-l)+lUn^) 

- m chx ^Cp (l - ln^ + R s ln(48) 
= m ch xRs(ln^-ln^ l -ln^-l), 


which, together with the thermal state equation of a perfect gas, 
results in 


dSi 


irr,cyl 


df 


= m chx Ks In 


T c p°V°\ 
pVcyl TO J 


= -rrichxRs In 


. m cyi 
m° 


+ 1 


where m cy i represents the gas mass content of the cylinder and 
m° represents a constant. Cycle integration, together with 
m chx = -dm^i/dt and 



(49) 


leads to 

V.cyi = j dS in . cy i = R s ^ln^ + 1 j dm cyl = 0. (50) 

The entropy production rate in the cold heat exchanger is 
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. m [hN Cp(l - In*) - m cte c p (If-In^) 

/ T" T" \ 

= rii chx c p ^l --^-+ln^J. 

(51) 

Cycle integration, together with Eqs. (10), (27), (29), and (34), as 
well as e = e e and e c = ioe, leads to 


_ CpPmax ( ( T 0 + l) g _ | 1 + £ \ 

R S T C V d + C e ^l-T 0£ ' 1 -T 0 e/ 


(52) 


The entropy production rate in the pulse tube is given by 

= (m ch x - ™ hhx ) (cv(ln^ - l) + K s ln^) 

+ m chx ^c p ^ - ln + Rsin jfiJ 
~ m hhx ^c p - ln^ + Rs\n^j 
= (m ch x - m hhx ) ^v^n^-l^)+Rs^n^+ln^jj 
+ -Wp ({g - ln^ - m hhx c p ^ , 

(53) 

which, together with m chx - m hhx = dm pt /df, the thermal state 
equation of a perfect gas, and R s = c p -c v , becomes 


_ c pPmin ^p r d(l — T 0 £ ') (, 1 + £ ( T 0 + l) fN \ 

’ hhx " RsT c r d (l — T 0 e) + 1 V 1 -T 0 e 1 + e / 

3 

+ Cp J m hhx - 1 - ln^j df. (56) 

2 1 1 

For the entropy production rate in the regenerator follows 


%* - + «n^) + ^(cp(i 

-ini) +R, |n |,), 

(57) 

in which all terms except of p and m hhx are constant. Thus, its cycle 
integration leads to 


^irr.reg ~ ( r ^^irr.cyl 


= ( c v ( ln^p- 1 ) +C V 


- In 


i re g 


+ R,ln^ 


)/ 


dmi 


hhx 


+ R,j m hhx ln^df ( = ) R s £ m hhx ln^df. 


JO 


(58) 


With Eq. (18) Eq. (58) becomes 
r d 


Si 


P (33) 

irr ' reg - (r c - l)(r d + 1) TVeg J *V dP = °- 


(59) 


dSj 


irr,pt 


dt 


= (m chx - mhhx)c P ln ^ ( 


m° 


, 7 pt \ . (T" t T" t 
+ m chx c p | ln^j -m hhx c p l^-ln^ 


where the gas mass content of the pulse tube is m pt and m° is a 
constant. Cycle integration, together with Eq. (49), leads to 


Sirr.pt = c p / (m chx - m hhx )in^dt + c p / [ m chx f ^ - in ^ 


pt 


JO 


( th th 

licit. 


The entropy production rate in the hot heat exchanger is 


jo 


(54) 


dSjrr.hhx _ ■ fT hhx | n ^hhx') ™ 

^ ~ m hhx c P ( y h3 J ™hhx c ] 


mhhxCp(-y^-l-ln-y 1X 


(55) 


5. Results 

A virtual pulse tube engine having the characteristics r c = r d = 1, 
Pmin^d = 1 J, and T c = 300 K has been simulated. It uses the perfect 
gas Helium with k = 1.66 as working fluid. The minimal pressure is 
p min = 1 bar and the displacement volume is i/ d = 10 cm 3 . The 
thermodynamic cycle was calculated for temperature ratios be¬ 
tween 1 (T h = 300 K) and 5 (T h = 1500 K). 



Cycle integration, together with Eqs. (10), (20), (22), and (34), as 
well as e = e e and e c = rof, leads to 


Fig. 6. Analytically calculated pV-work and amounts of heat transferred per cycle in 
the cylinder and the pulse tube as functions of the temperature ratio. 
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Fig. 7. Analytically calculated amounts of heat transferred per cycle in the cold and hot 
heat exchanger as functions of the temperature ratio. 


5.2. Calculated pV—work, transferred heat, thermal efficiency, and 
entropy production 

Figs. 6 and 7 show the pV—work, calculated with Eq. (39), and the 
amounts of heat transferred per cycle in the engine components, 
calculated with Eq. (46). All energies are zero at t = to = 2. For t < to, 
the pV—work is positive and work has to be delivered to the engine, 
which is thus not acting as a prime mover. In this case heat is added 
to the working gas in the cold heat exchanger and rejected from it in 
the hot heat exchanger. The engine operates as a heat pump. For 
t > to, the pV—work is negative and the engine operates as a prime 
mover. Heat is added to the working gas in the hot heat exchanger 
and rejected from it in the cold heat exchanger. The amount of heat 
transferred per cycle in the pulse tube is much smaller than in the 
other components. According to Fig. 6, the heat transferred per cycle 
in the pulse tube has an extreme value at t = to. Thus, 

«PminV , d(r d (to - 1) - 1) ' a rcn , 

(«-l)(r fl + l).„ -° (60) 

With Eq. (60) the minimal temperature ratio above which the en¬ 
gine provides a work output results in 


dQpt 

dr 


to 


td + 1 
td 


The second derivative of Qp t in t at t = tq is 


(61) 


d 2 Qpt 

dr 2 


T=T 0 


KPmmVqid 4r 2 (r c -l) + 2r c r d (r c + 2) + r 2 

r c (K—1) (r d + l) 2 (r d (2r d + 3) + l)(2r d (r c + l) + r c )' 

(62) 


The sign of Eq. (62) is governed by 
a = 4r d (r c - 1) + 2 r c r d (r c + 2) + r 2 . 
If r c < 1 and 


r 2 + 4r d + 2 - (r d + 1) 
2r d + 1 


r c < 2r d 


a is smaller than zero and thus d 2 Qp t /dT 2 | r=Io < 0. In this case, at 
t = to, a local maximum of Qp t is located and integrated over one 
cycle heat is removed from the working gas in the pulse tube. 

If r c < 1 and 


Jrj + 4r d + 2- (r d + l) 
r ‘ - 2r -- --■ 

a is equal to zero and thus d 2 Qpt/dT 2 | I=To = 0. In this case, inte¬ 
grated over the cycle, no heat is exchanged between the working 
gas and the pulse tube material. In all other cases, integrated over 
the cycle, heat is added to the working gas in the pulse tube. This 
situation allows for an addition of heat to the working gas at 
temperatures lower than Th, enabling the engine to convert heat 
energy from sources with decreasing temperature. Only this case is 
of practical interest. 

The thermal efficiency is calculated with 


Vth = 


W pV 

Qin 


(63) 


where Qi n represents the heat input. According to Figs. 6 and 7, for 
t = t 0 , heat is added to the working gas in the cylinder at the 
temperature T c , in the hot heat exchanger at the temperature Th, 
and for a > 0 in the pulse tube at temperatures between T c and Th. 
The heat added to the working gas in the cylinder, in general, is 
provided from the surroundings and does not have to be provided 
by the heat source. Taking into account only the heat delivered by 
the hot heat exchanger at the temperature Th, the heat input is 
calculated with Qj n = Qhtix- Considering in the case of a > 0 also the 
heat added to the working gas in the pulse tube at temperatures 
between T c and Th, the heat input is calculated with Qj n = Qhhx + Qpt- 
In both cases, for t — > tq, Eq. (63) becomes finite, equaling to 


k— 1 

Vth = x(r d + 1) ’ 


(64) 


which is only possible if to > 1. 

In Fig. 8, the thermal efficiency of the pulse tube engine, 
calculated with Eq. (63), is shown for temperature ratios between 
t = to and t = 5. At t = to, the efficiency is highest, given by Eq. (64), 
and decreases monotonously for t > tq. The efficiency calculated 



Fig. 8. Analytically calculated thermal efficiency of the pulse tube engine as a function 
of the temperature ratio. 
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Fig. 9. Analytically calculated entropy production per cycle in the cold and hot heat 
exchanger as well as in the pulse tube as functions of the temperature ratio. 

with Qin = Qhhx is higher than the efficiency calculated with 
Qin = Qhhx + Opt, but in the case of Qj n = Qhhx + Opt. heat is added to 
the working gas at temperatures between T c and Th in the pulse 
tube. 

The decreasing nature of the thermal efficiency for x > i 0 is 
based on the working principle of the pulse tube engine. This be¬ 
comes obvious regarding Fig. 9, which shows the entropy produc¬ 
tion in the cold heat exchanger, the pulse tube, and the hot heat 
exchanger, calculated with Eqs. (52), (54), and (56). For x + to, 
entropy is produced in all components, which is caused by mixing 
of gas at different temperatures. Since a pV —work output is only 
provided for x > to, the pulse tube engine is intrinsically irrevers¬ 
ible. This was also concluded from experimental investigations of 
the pulse tube engine by Yoshida et al. [13], 

5.2. Influence of operating and design parameters on the engine 
cycle 

The time-dependent pressure of the engine during compression 
and expansion can be calculated with Eqs. (21) and (28). The results 
are used to analyze the engine cycle. 



Fig. 10. Analytically calculated pV-diagram of the virtual pulse tube engine for r c = 1, 
r d = 2 and temperature ratios of 1, 2, and 3. 


The influence of the temperature ratio on the thermodynamic 
cycle was investigated. Fig. 10 shows the resulting pV— diagram of 
the virtual pulse tube engine for r c = 1, rj = 2 and temperature 
ratios of 1, 2, and 3. The minimal temperature ratio for engine 
operation is to = 1.5 (Th = 450 K). The slopes of the isothermal 
compression and expansion increase with an increase of t. For t = 1, 
the expansion stroke is on a lower pressure level than the 
compression stroke. In this case, pV— work has to be delivered to the 
engine. For x = 2 and x = 3, the expansion stroke is on a higher 
pressure level than the compression stroke and the engine provides 
a pV— work output. The intensity of the isochoric pressure rise at 
TDC is highly dependent on x. In the case of x = 1, the gas in 
the pulse tube is cooled at the end of the compression stroke, since 
its compression end temperature is higher than the pulse tube’s 
wall material temperature. For higher temperature ratios, the 
compression end temperature is lower than the pulse tube’s wall 
material temperature resulting in a heat addition to the working 
gas at TDC. The same happens at BDC, but with an opposite sign of 
the heat flux in the pulse tube. If the temperature ratio is too small, 
the working gas expands below the pulse tube’s wall material 
temperature and takes on heat at BDC. Having higher temperature 
ratios, the gas temperature in the pulse tube does not fall under the 
pulse tube’s wall material temperature and heat is rejected at BDC. 

Next, the influence of the compression ratio on the engine cycle 
is analyzed. In Fig. 11, the pV— diagram of the virtual pulse tube 
engine for ra = 2, x = 2 and compression ratios of 0.001,1, and 2 is 
shown. As expected from the theory, the pressure rise during 
compression increases for higher values of r c . For all compression 
ratios, the isothermal expansion is on a higher pressure level than 
the isothermal compression. Thus, the engine provides a work 
output for all values of r c > 0. The surface area enclosed by the pV— 
curve is comparable for all compression ratios. For r c -> 0, the 
temperature change during compression and expansion becomes 
zero. According to Eq. (54), the entropy production in the pulse tube 
becomes singular, and thus, the energy dissipation becomes infin¬ 
ite. This hampers an operation of the engine for r c —> 0. The 
calculated work output for r c —» 0 is only a theoretical possible 
value for Qj n —► °°. 

According to Eq. (39), the pV —work is only of minor order 
dependent on r c . To show this, the pV— work is plotted in Fig. 12 as a 
function of the compression ratio, for r^ = 2 and temperature ratios 
of 2, 3, 4, and 5. The values are normalized to the value for r c = 0. 
The dependency of the pV— work on r c is evident. For small 



Fig. 11. Analytically calculated pV-diagram of the virtual pulse tube engine for r = 2, 
r d = 2 and compression ratios of 0.001,1, and 2. 
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Fig. 12. Analytically calculated normalized pV -work as a function of the compression 
ratio and entropy production in the pulse tube as a function of the compression ratio of 
the virtual pulse tube engine having a distribution ratio of 2 and temperature ratios of 
2, 3, 4, and 5. 



Fig. 14. Analytically calculated pV-diagram of the virtual pulse tube engine for t = 2, 
r c = 1 and distribution ratios of 0.5,1.2, and 3. 


temperature ratios, the influence of r c is insignificant. Nevertheless, 
an increase of r c slightly reduces the pV —work of the ideal pulse 
tube engine. The influence of the compression ratio on the pV— 
work increases with the temperature ratio. On the other hand, 
according to Eq. (54), an increase of r c decreases the entropy pro¬ 
duction in the pulse tube. This is also shown in Fig. 12. Thus, the 
power output of the real pulse tube engine should have an optimal 
compression ratio at which the power output is best. This peak 
value should be located in the region of high dSj rriPt /dT around 
r c = 1, which was proven numerically in Ref. [23], There, the used 
compression ratio r* is defined as r* = r c + 1. 

The thermal efficiency, calculated with Eq. (63), is shown in 
Fig. 13. With increasing compression ratios up to r c = 1, the thermal 
efficiency increases rapidly and remains constant for higher 
compression ratios. This strengthens the result of a peak pV —work 
around r c = 1. 

Finally, the distribution ratio is studied with respect to its in¬ 
fluence on the engine cycle. Fig. 14 shows the pV— diagram of the 
virtual pulse tube engine for r c = 1 , t = 2 and distribution ratios of 
0.5, 1.2, and 3. For ra = 0.5, the engine does not provide a work 
output. Increasing the distribution ratio enables the cycle to 


convert heat into work. An increase of the distribution ratio in¬ 
creases the regenerator volume and decreases the pulse tube vol¬ 
ume. This directly influences the thermodynamic asymmetry 
caused by the arrangement of pulse tube and regenerator. A higher 
value of rd increases the mass flux passing the pulse tube and the 
hot heat exchanger, but decreases the mass content experiencing a 
thermal relaxation in the pulse tube. Thus, an optimal value of r d , 
maximizing the pV— work, exists. The linear component of the first 
derivative of the Eq. (39) in ra is 


dW pV j 

dr d 


Pmin^d— 


-P-(T+P 

( r d + l) 3 


which becomes zero for 


(65) 


^ + 1 f CC\ 

r d,opt ~ 'i — \ ' (66) 

For different temperature ratios, the pV —work is calculated and 
shown in Fig. 15. The strong dependency of the pV —work on the 
distribution ratio is evident. The influence of the regenerator on the 
cycle becomes obvious. The lower the available temperature ratio, 



Fig. 13. Analytically calculated thermal efficiency as a function of the compression 
ratio of the virtual pulse tube engine corresponding to the conditions of Fig. 12. 



Distribution ratio 

Fig. 15. Analytically calculated pV-work as a function of the distribution ratio of the 
virtual pulse tube engine having a compression ratio of 1 and temperature ratios of 1.5, 
2, 2.5, 3, and 3.5. 





































296 


S. Moldenhauer / Energy 62 (2013) 285-299 



Fig. 16. Analytically calculated thermal efficiency as a function of the distribution ratio 
of the virtual pulse tube engine corresponding to the conditions of Fig. 15. 

the higher the value of ra and thus, the higher the necessary 
regenerator gas volume to enable a pV —work output. As expected, 
exceeding the optimal distribution ratio leads to a decrease of the 
pV —work output. A similar result was obtained numerically in Ref. 
[23]. Since according to Eq. (66) the peak value is a function of the 
temperature ratio, a careful design of the pulse tube engine with 
respect to the distribution ratio is necessary. The found relation in 
Eq. (66) delivers an approximation of the optimal distribution ratio. 
In contrast, the thermal efficiency, calculated with Eq. (63) and 
shown in Fig. 16, does not have a peak value at any distribution 
ratio. To increase the thermal efficiency, the distribution ratio has to 
be small. According to Eq. (61), the minimal temperature ratio tq 
above which the engine provides a work output is a function of ra. 
Thus, the lowest possible value of rd is limited by the available 
temperature ratio. 

5.3. Comparison of the theoretical results with experimental data 

In Refs. [4,24] a working pulse tube engine was presented. The 
experimental test setup is used to prove the predictions of the 
analytical model. 



Temperature ratio 


Fig. 17. Fits of the measured net power outputs as functions of the temperature ratio 
for filling pressures of 6 bar, 8 bar, and 10 bar, an operating frequency of 8 Hz and a 
cooling temperature of 25 °C. 


The gross power P g of the pulse tube engine is related to the pV— 
work through P g = fW p y. The net power output of the engine is 
P n = Pg-Pi i n t, where the total internal power loss Pi j nt , in general, is 
a function of at least/, pf, T, geometry parameters such as r c , rd, and 
I'd as well as working gas properties. The geometry of the engine 
(r c = 0.38, ra = 1.71, and Vd = 18 cm 3 ) is fixed. Thus, the internal 
power loss is reduced to being a function of f p f, T, and the working 
gas properties. The main influence of the temperature on the in¬ 
ternal power loss is due to the temperature dependency of the 
working gas properties. Assuming a perfect gas and neglecting the 
remaining secondary dependency of the internal power loss on the 
temperature, the internal power loss is only a function of the fre¬ 
quency / and the filling pressure pf. 

The net power can thus be fitted with 

Pn = FfW P V ~ Pl,int, (67) 

where p represents a constant fitting parameter. 

For a fixed frequency and filling pressure, the factor p and the 
total internal power loss Fl int can be obtained from the experiment. 
To do so, the net power output was measured for filling pressures of 
6 bar, 8 bar, and 10 bar and frequencies from 4 Hz to 12 Hz. For each 
couple of (/,Pf), the engine was motored at heat input temperatures 
between 40 °C and 600 °C. The minimum pressure is considered to 
be independent of the temperature distribution. In reality, due to 
the thermal relaxation, the temperature distribution slightly in¬ 
fluences the minimum pressure and the pressure amplitude. The 
experimentally observed minimum pressures are p m jn = 5.15 bar 
for a filling pressure of 6 bar, p m j n = 7.00 bar for a filling pressure of 
8 bar, and p m i n = 8.55 bar for a filling pressure of 10 bar. 

The measured net power P n is plotted versus the temperature 
ratio i = 7h/T c and adapted with Eq. (67). The fitting parameters are 
P and Fl int. For f = 8 Hz, the result is shown in Fig. 17. 

Almost independent of the filling pressure, the net power be¬ 
comes zero for x = 1.62 ± 0.05. The respective minimal heat 
input temperature above which the engine is working is 
Th = 210 °C ± 15 K. The minimal temperature ratio above which the 
engine provides a gross power output was found with Eq. (61) to be 
i 0 = (id + 1 )/i'd. The distribution ratio of the engine is r<j = 1.71 and 
thus ro = 1.59. This value is within the range of the experimentally 
observed value and differs from the measured mean value by about 
2%. The result is highly interesting, since simplifications such as 



Frequency [Hz] 

Fig. 18. Internal power losses as functions of the operating frequency for filling 
pressures of 6 bar, 8 bar, and 10 bar, obtained from the fit of the measured net power in 
Fig. 17 with Eq. (67). The cooling temperature is 25 °C. 
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Fig. 19. Comparison of the analytically calculated and the measured minimal tem¬ 
perature ratio above which the engine provides a power output as functions of the 
distribution ratio. The experiments are done for filling pressures of 6 bar and 10 bar as 
well as operating frequencies of 5 Hz and 7 Hz. The cooling temperature is 25 °C. 

isothermal compression and expansion were made for the analyt¬ 
ical approach. It seems that the property of the pulse tube engine to 
have the calculated and measured minimal temperature ratio to, is 
a fundamental characteristic of the engine. The analytical model is 
robust enough to detect this property. The model predicts that the 
minimal temperature ratio is only a function of the distribution 
ratio. In the experiment, a weak dependency of to on the filling 
pressure was observed. 

The fitting factor f is obtained to be /- = 0.159 ± 0.016. The in¬ 
ternal power loss, also obtained from the fitting procedure, is 
plotted in Fig. 18. Up to an operating frequency of about 8 Hz, the 
internal power loss is very low (< 0.5 W). After exceeding this 
frequency it increases significantly. Fitting the internal power loss 
with the function f\j n t = Q / delivers an exponent of l ~ 5. Me¬ 
chanical power loss due to deformation of the regenerator mesh is 
proportional to/and flow friction power loss is proportional to f. 
Thus, the observed internal power loss does not have its origin in 
mechanical or flow friction. The thermodynamic situation in the 
engine seems to change for higher frequencies, leading to the 
observed strong increase of the internal power loss. 



Temperature ratio 

Fig. 20. Comparison of the measured thermal efficiency for operating frequencies 
between 4 Hz and 12 Hz with the analytically and numerically calculated thermal 
efficiency as functions of the temperature ratio. The measured efficiencies are available 
up to a heat input temperature of 600 °C. The filling pressure is 10 bar and the cooling 
temperature is 25 °C. The respective Carnot cycle efficiency is also shown. 


In order to further investigate the validity of Eq. (61 ), the size of 
the regenerator of the working pulse tube engine of [4,24] has been 
modified. The regenerator is a stack of wire mesh screens having a 
mesh number of 50, a wire diameter of 200 pm, and a screen 
diameter of 20 mm. The length of the regenerator was varied in 
steps of 10 mm between 50.5 mm and 110.5 mm resulting in dis¬ 
tribution ratios between 0.64 and 1.33 for the 50.5 mm and 
110.5 mm regenerator lengths, respectively. For each regenerator 
size, the minimal temperature ratio above which the engine pro¬ 
vides a net power output was measured at filling pressures of 6 bar 
and 10 bar for operating frequencies of 5 Hz and 7 Hz. According to 
Fig. 18, the internal power loss at frequencies lower than 8 Hz is 
negligible, so that P g = P n . For all the distribution ratios investi¬ 
gated, the minimal temperature ratio above which the engine 
provides a gross power output was calculated with Eq. (61). The 
results are shown in Fig. 19. The relative error of the measured 
cooling temperature is estimated to be 1%. Due to the difficulty in 
the adjustment of the correct heater temperature, its relative error 
is about 10%. Thus, the error of the measured minimal temperature 
ratio is about 11%. The error bars are also shown in Fig. 19. 

For all filling pressures and frequencies, the analytically calcu¬ 
lated minimal temperature ratios are within the range of uncertainty 
of the measured values. The calculated and measured trends also 
match. Eq. (61) does not predict a dependency of to on the filling 
pressure and the operating frequency, which is not in accordance 
with the observations in the experiments. As shown in Fig. 19, for a 
fixed frequency, a higher filling pressure leads to a higher minimal 
temperature ratio, and for a fixed filling pressure, a higher frequency 
leads to a lower minimal temperature ratio. For the working engine, 
the decrease of to with an increase in the operating frequency is a 
result of an increase in the phase delay between pressure swing and 
heat flux in the pulse tube (thermal lag) at higher frequencies. A 
higher filling pressure on the other hand increases the heat flux and 
thus decreases the thermal lag. This leads to an increase in to at 
higher filling pressures. 

In Fig. 20, the measured thermal efficiency of the pulse tube 
engine is compared with the analytically calculated thermal effi¬ 
ciency according to Eq. (63) in combination with Eqs. (39) and (46) 
and with the respective Carnot cycle efficiency. The analytically 
calculated thermal efficiency is an upper limit of the efficiency of 
the pulse tube engine. It is highest for t = t 0 and decreases with an 
increase of t. The measured thermal efficiency of the working 
pulse tube engine does not show this behavior. It is zero for t = to 
and increases with an increase of t. Since the gross efficiency 
cannot exceed the analytically calculated (ideal) thermal effi¬ 
ciency, it has to decrease above a certain temperature ratio T opt , 
which was too high for an experimental investigation in this work. 
To show the existence of T op t, the numerical simulation model of 
the pulse tube engine presented in Ref. [24] was used. The pulse 
tube engine was simulated for heat input temperatures up to 
1000 °C (t = 4.27) and the thermal efficiency is calculated 
numerically. The result is also shown in Fig. 20. The numerically 
calculated thermal efficiency reaches its peak value at about 
T op t ~ 3. For higher temperature ratios, the thermal efficiency 
decreases slowly. Further numerical calculations for higher tem¬ 
perature ratios show that the numerically calculated thermal 
efficiency never reaches the value of the analytically calculated 
thermal efficiency. For a heat input temperature of 1500 °C 
(t = 5.95), the numerically calculated thermal efficiency is 5.5%, 
which is lower than the analytically calculated thermal efficiency 
of 5.8%. At t = to, the numerically calculated thermal efficiency is 
2%. This indicates the existence of a finite efficiency at t = to as 
predicted by the analytical model, but the power output of the 
working engine for t -» to is too small to allow an experimental 
determination of the respective efficiency. 
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For the working engine, the thermal efficiency increases with 
higher operating frequencies. The reason is the increase in the 
phase delay between pressure swing and heat flux in the pulse tube 
(thermal lag). For higher frequencies, the time at which heat is 
added to the working gas and removed from the working gas in the 
pulse tube gets closer to TDC and BDC, respectively. The amount of 
heat transferred due to thermal relaxation in the pulse tube in¬ 
creases and the cycle gets closer to the theoretical cycle. 

Compared to the Carnot cycle efficiency, the efficiency of the 
pulse tube engine is small. The highest theoretical value of the 
pulse tube engine’s efficiency is for x = to with a value of 14.7% only 
40% of the respective Carnot cycle efficiency. Thus, the irreversible 
nature of the pulse tube engine is clearly evident. Regarding the 
highest measured efficiency of the working pulse tube engine, its 
irreversible operation becomes even more obvious. The experi¬ 
mentally achieved efficiency of about 8% at i = 2.9 amounts to only 
12% of the respective Carnot cycle efficiency. 

6. Conclusions 

The present work represents an analytical approach on the basic 
characteristics of the pulse tube engine. The main results are 
summarized as follows: 

1. The pV— work of the pulse tube engine can be estimated with 
Eq. (39). 

2. According to Figs. 6 and 9, entropy is always produced if heat is 
converted into pV —work and vice versa. Thus, the irreversible 
nature of the pulse tube engine (i > to) and the pulse tube heat 
pump (t < to) is proven theoretically. 

3. The ideal pulse tube engine starts to operate at a temperature 
ratio of t > to > 1, where according to Eq. (61 ), to = (ra + 1 )/r c | is 
determined by the distribution ratio. The pV— work output is 
zero for t = to and increases monotonously for t > to, whereas 
the thermal efficiency has its maximum at t = to and decreases 
monotonously for t > to. This is a result of the strong increase of 
the entropy production and only a moderately increasing work 
output. Thus, the pulse tube engine is preferably suited for low 
temperature ratios t, which are bound to the condition t > to- 

4. The theoretically found minimum temperature ratio to above 
which the pulse tube engine is working is confirmed experi¬ 
mentally. The measured value amounts to t 0 = 1.62 
(Th = 210 °C), with an absolute error of At = 0.05 (ATh = 15 K). 
The value calculated with Eq. (61) amounts to to=1.59 
(Th = 210 °C) and is thus within the error range of the measured 
value. Since Eq. (61) is obtained from an idealized analysis of 
the pulse tube engine and the measurement delivers the same 
result, the minimal temperature ratio to seems to be a basic 
property of the pulse tube engine and the developed analytical 
model is robust enough to determine this basic engine prop¬ 
erty. Moreover, for different values of r { j, the analytically 
calculated trend of to is confirmed experimentally. In contrast 
to the predictions of Eq. (61), the experimental investigations 
disclose that to increases for higher filling pressures and de¬ 
creases for higher operating frequencies. 

5. For the ideal engine, the heat added to the working gas in the 
cylinder is converted into pV— work. The entire heat added to 
the working gas in the hot heat exchanger and the pulse tube is 
rejected in the cold heat exchanger. This effect enables a pulse 
tube engine with a non-isothermal cylinder to provide a cool¬ 
ing power in the cylinder. The pulse tube engine can thus act as 
a thermal driven refrigerator, absorbing heat in the cylinder. 
Even if the thermal energy applied in the hot heat exchanger 
and the pulse tube is not converted into pV— work, it is essential 
for the engine to operate. 


6. Since heat can be added to the working gas in the pulse tube at 
temperatures between T c and Tj,, the pulse tube engine is able 
to use heat sources at different temperatures or absorb heat 
from a fluid flowing along the pulse tube from its hot to its cold 
side. It is theoretically possible to design a pulse tube engine, 
which cools down a fluid, beginning at the hot heat exchanger 
at the temperature Th, along the pulse tube to the cold heat 
exchanger at the temperature T c . This allows for the removal of 
a larger amount of heat from a thermal energy carrying fluid 
compared to e.g. the Stirling engine, which due to its regen¬ 
erative cycle is only able to use heat at the temperature of the 
hot heat exchanger Th. 

7. The pV— work and the efficiency of the pulse tube engine are 
strongly dependent on the distribution ratio. The pV —work 
output is highest for a certain value of r^, which can be esti¬ 
mated with Eq. (66) to rd.opt ~ (t + 1)/(t— 1). The thermal ef¬ 
ficiency decreases monotonously with increasing distribution 
ratios. 

8. An upper limit for the efficiency of the pulse tube engine is 
found in Eq. (63) and Eq. (64). The measured and numerically 
calculated gross efficiency never exceeds this limit, which is 
highest for t = to, and decreases monotonously for t > to. The 
gross efficiency of the working engine is zero for t = to. It in¬ 
creases up to a certain temperature ratio T opt , which is found to 
be about T opt = 3, and decreases for higher temperature ratios. 
Due to the irreversible nature of the pulse tube engine, the 
efficiency of the investigated ideal engine never exceeds 40% of 
the respective Carnot cycle efficiency. The experimentally 
achieved highest efficiency [24] is only about 12% of the 
respective Carnot cycle efficiency. 

The listed properties of the pulse tube engine give a pathway for 
its potential applications and engineering challenges. 
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